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1. INTRODUCTION 

The use of vertical-incidence Doppler velocity in 
addition to radar reflectivity may yield information 
on drop size distribution and therefore result in 
better rainrate estimates. Doppler velocity can 
provide useful information on the raindrop size 
distribution. Doppler velocities from a zenith- 
pointing radar represent the sum of the mean 
reflectivity-weighted hydrometeor fallspeed and the 
vertical air motion. Dual-parameter rain estimation 
methods using the Doppler velocity, require that the 
latter can be removed, or is negligible. Atlas et al. 
(1972) derived relations between Doppler velocity, 
reflectivity, and rain rate assuming an exponential 
size distribution for rain. Ulbrich (1994) expanded 
on this work by deriving the relation between the 
Doppler velocity and the reflectivity assuming a 
Gamma size distribution. This distribution provides 
a more realistic representation of the small rain 
drops. To get accurate information on raindrop size 
distributions with the above method, the air motions 
must be removed from the observed Doppler 
velocities. 

The purpose of this paper is a preliminary 
examination of raindrop size distributions using 
dual-parameter (reflectivity, Doppler velocity) 
airborne radar measurements from nadir-pointing 
X-band ER-2 Doppler radar (EDOP) located on the 
high-altitude ER-2 aircraft (Heymsfield et al. 1996). 
EDOP provides similar measurements to ground- 
based vertically pointing radars except that 
extensive regions can be covered in short times; 
ground-based vertical-incidence radars provide 
time-height measurements which and in all but a 
few precipitation systems time is a poor substitute 
for distance. The ER-2 participated in three recent 
Tropical Rain Measuring Mission (TRMM) field 
campaigns: Texas and Florida Underflights 
Experiment (TEFLUN-A and TEFLUN-B), and 
TRMM-Brazil (Amazonia). TEFLUN was conducted 
for the purpose of validation of the rain 
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measurements from the TRMM satellite. In addition, 
flights were made during the Convection and 
Moisture Experiment (CAMEX-3) which was aimed 
at hurricane studies. The current data sets from 
EDOP provide the opportunity to examine the 
general properties of the vertical structure of the 
raindrop size distributions in these different 
precipitation regimes. 

EDOP samples Doppler velocity, v D , and reflectivity 
factor, Z, with a gate spacing of 37.5 m (vertical 
resolution) and at 0.5 s intervals (100 m along-track 
sampling frequency). The antenna beamwidth is 
about 3° which provides a spot size of about 1 .2 km 
at the surface. Calibration accuracy of EDOP is 
about 1-2 dB. 

2. APPROACH 

The radar reflectivity factor 2 in mm 6 nrr 3 can be 
given by (e.g. Ulbrich 1992): 

(1) 

where D (cm) is particle equivalent spherical 
diameter, N(D) in m^cm* 1 is number of raindrops 
per unit volume per unit size interval, D 0 is the 
median particle diameter, the Gamma size 

distribution is assumed N(D) = N 0 D^e~ AD . 
where p is the shape parameter, and a =A D 0 = 3.67 
+ p + ia°- 3( 9 ^ } . A/ 0 ,( p>= 6 A x 10* exp(3.2p)) has 
been shown observationally and theoretically by 
Ulbrich (1983). The integration limits D m , n and D^ 
are assumed 0-»°° , though several papers have 
dealt with truncated limits. The reflectivity-weighted 
Doppler velocity v z (ms 1 ) is given by: 

Jn m “ v(D)N(D)D 6 dD 

v = mln (2) 

Z f % ma N(D)D 6 dD 

where v(D) is the drop fallspeed distribution. 
Furthermore, Ulbrich (1994) has derived the 
reflectivity-weighted Doppler velocity v z as a 
function of p by assuming the Gunn-Kinzer 
fallspeed relation at ground level (1013 mb), i.e.. 
v(D) = 9.65 - 10.3 e' €D . Assuming that vertical air 
motions are neglible, Eq. 2 by elimination of 
parameters can be given as a function of u and 2: 
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v z (p) = 9.65 




-10.3 1 + 6 


N 0 (ii)\o b rn+v) 


(3) 

This relation applies to ground level (1013 mb); v z 
at other altitudes is obtained by multiplying it by 
[pc/p] 044 where p and p 0 are the air density at the 
surface and measurement height, respectively. Eq. 
3 has a major advantage over power law relations 
used extensively in the literature: v z is asymptotic 
with increasing reflectivity rather than increasing 
monotonically with Z The rain rate R can also be 
computed as a function of p (Ulbrich 1992): 

/?(/J) = F(p)Zv- 3 - 48 , (4) 

where F is an analytic function of p. 


Figure 1 shows the dependence of v z on fusing 
Eq. 3; a commonly used power law relation is also 
plotted. At 40 dBZ, the variation of p from -2 to 6. 
results a fairly substantial 4 ms* 1 variation in v z .. 
Figure 2 shows the Gamma size distribution 
corresponding to v z = 6, 8 ms* 1 and p =-2 to 4. The 
resulting range of generated curves are typical of 
reflectivities observed in EDOP data sets. Larger v z 
for a given p implies larger D 0 and a larger Z 



Figure 1. v z - Z relation as a function of p. 


The approach here is to fit Eq. 3 to the EDOP 
vertically pointing Doppler velocity v D and Z 
observations to obtain pin the Gamma distribution. 
The main assumption in applying Eq. 3 to EDOP 
observations is that stratiform regions are fairly 
extensive and that the mean vertical air motions 
over this region are negligible. Eq. 3 is modified 
with an additive constant because widespread 
vertical motions can be present. Thus, if one has 
large sample of v D and Z measurements from the 
same region, then p can be obtained by fitting the 
altitude-adjusted observed v D with the following: 

z \ 0.44 

— v D =C+v z {n), (4) 

\Po J 

where v z is obtained from Eq. 3 and C is a constant 
which consists of the mean mesoscale vertical 
motion plus biases caused by incorrect removal of 
aircraft motion. This fitting method provides a 
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Figure 2. Gamma distribution for two values of v z , 
assuming various values of the shape parameter p. 


global estimate of p and C within some defined 
height interval, not a point-by-point estimate of the 
values. Besides the inhomogeniety of the size 
distribution over the analyzed region, a major 
source of error in this approach is the presence of 
small scale vertical air motions. In addition, small 
scale aircraft motions are not always fully removed 
resulting in an error in the v D used in the 
calculation. 


3. EXAMPLE. 


We provide an example of Eq. 4 applied to data 
(Fig. 3) from an ER-2 flight on 13 January 1995 
covering a stratiform region trailing a squall line off 




the Gulf coast of the U.S. The rain region between 
the surface and 400 m below the freezing level, is 
divided into 525 m intervals and then fitted using a 
non-linear curve fitting routine. The vertical axis 
corresponds to observed Doppler velocities 
adjusted to the surface, with a range of values from 
about 3 to 9 ms* 1 . The reflectivities in the data set 
range from about 20 to 45 dBZ; no attenuation 
correction has been applied to this flight line. About 
10,000 points went into the fit and while there is 
scatter of the points, the general shape of the 
cluster of points is fitted reasonably well by the 
curves. The fitted curves in Fig. 3 indicate that 
p ranges from -0.84 at the top layer (1.5-2 km) to 
1.2 at the bottom layer (0.4-0. 9 km), implying a 
nearly exponential distribution. The calculated C 
ranges from about +0.6 ms* 1 at the top to +1.3 ms* 1 
in the bottom layer which suggests either a 
mesoscale downdraft or aircraft motions. For the 
same v z , a larger p implies a shift toward larger 
droplets. This would be consistent with 
subsidence-induced evaporation as was present 
with this case. The downward increase in C is 
consistent with the subsiding rear inflow observed 
with the squall line; this is probably not caused by 
aircraft motions which would affect all levels 
uniformly. 
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Bottom Top Npts /i C N q 

(km) (km) (ms -1 ) (m~ 3 cm~ 1 ~'‘) 

0.30 0.83 10408 1.2 1.3 3.4E+06 

0.90 1.43 10408 0.3 1.0 1.5E+05 

1.50 2.03 10408 -0.8 0.6 4.3E+03 


Figure 3. Observed measurements from squall line 
stratiform on 13 January 1995 fitted with Eq. 4. 


5. CONCLUSIONS. 

This paper has examined an approach for 
estimating parameters of Gamma rain size 
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distribution using only Doppler velocity and 
reflectivity measurements from an airborne nadir- 
pointing Doppler radar. The approach is best 
suited for widespread uniform stratiform regions 
without significant small-scale vertical motions 
present. Although this approach has some 
limitations when compared with other 
multiparameter approaches (e.g., dual-wavelength, 
dual-polarization, Doppler spectrum), it provides 
useful information for understanding general 
characteristics of the drop size distributions, rain 
rates, and liquid water content. Additional cases 
are being examined with the above approach. In 
particular, this method may prove useful in 
estimating drop size distributions and rain rates in 
feeder bands of hurricanes sampled by EDOP. In 
addition, subtracting the best estimate of v z (p) from 
v D can provide an improved estimate of local air 
motions. 
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